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ABSTRACT 

Milliarcsecond VLBI maps of regions containing 6.7 GHz methanol maser emission have lead to the 
recent discovery of ring-like distributions of maser spots and the plausible hypothesis that they may 
be tracing circumstellar disks around forming high mass stars. We aimed to test this hypothesis by 
imaging these regions in the near and mid-infrared at high spatial resolution and compare the observed 
emission to the expected infrared morphologies as inferred from the geometries of the maser rings. 
In the near infrared we used the Gemini North adaptive optics system of Altair/NIRI, while in the 
mid-infrared we used the combination of the Gemini South instrument T-ReCS and super-resolution 
techniques. Resultant images had a resolution of ^150 mas in both the near-infrared and mid-infrared. 
We discuss the expected distribution of circumstellar material around young and massive accreting 
(proto)stars and what infrared emission geometries would be expected for the different maser ring 
orientations under the assumption that the masers are coming from within circumstellar disks. Based 
upon the observed infrared emission geometries for the four targets in our sample and the results of 
SED modeling of the massive young stellar objects associated with the maser rings, we do not find 
compelling evidence in support of the hypothesis that methanol masers rings reside in circumstellar 
disks. 

Subject headings: circumstellar matter - infrared: stars - masers - stars: formation 
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1. INTRODUCTION 

The formation of OB stars is an important challenge to 
modern astrophysics as they are responsible for many of 
the energetic phenomena in galaxies. But their large dis- 
tances, heavy obscuration and rapidity of evolution make 
studies of massive star formation difficult. One property 
that is unique to massive young stellar objects is the pres- 
ence of methanol masers in their e arly formative phases 
(|MenterilT99"il iCaswell et all 119951) . As such, they may 
hold clues to understanding the fundamental differences 
between low and high mass star formation. High angu- 
lar resolution observations of the masers allow us to map 
the neutral gas at scales of a few 10 AUs in the vicinity 
of massive young stars. However, because of their rela- 
tively sparse sampling, it is still unclear where and how 
the masers are formed. Two competing hypotheses have 
arisen; one stating that methanol masers are embedded 
in circumstellar tori or accretion disks around the mas- 
sive protostars (e.g., Norris et al. 1993), the ot her statin; 
that they may generally be tracing outflows (|De Buizc 
12001 . 

However, milliarcsecond scale very long baseline inter- 
fcromctry (VLBI) observations show a wide range of mor- 
phologies for the maser spot distributions at 6.7 GHz. 
They form simple structures (a single maser spot), or 
group into spots that lie in lines or arcs , or are dis- 
tributed randomly without any regularity (IMinier et al 
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19981 ). Linear structures of maser emission accompanied 



by clear velocity gradient s have led to the belief that they 
may be in edge-on disks (|Minier et al . 2000; Wals h et al.l 
19981). Other scenarios where outflows or shocks are gen- 
erating th e linear or arc-like d istributions have also been 
proposed (jDodson et al.ll2004f ). 

We have recently completed a survey of 31 sources 
at 6.7 GHz using Eur opean VLBI Network (EVN) 
(Bartkiewicz et al. 2009). In addition to the curved and 
complex morphologies observed in other samples, we 
have discovered for the first time nine sources (29 % of the 
sample) with ring-like maser distributions with typical 
sizes having major axes of 0'.'2-0'.'3. Though not apparent 
in the radio data, these ring-like structures strongly sug- 
gest the existence of a central stellar object, and prompt 
the obvious question: Are these maser rings tracing disks 
around massive protostars? 

Inspection of mid-infrared (MIR) data from the Spitzer 
IRAC maps, GLIMPSE and MIPSGAL 1 , revealed that 
all of our sources with ring-like morphologies coin- 
cide with unresolved MIR sources within one pixel 
in a GLIMPSE map (l'/2). Many sources also have 
bright near-infrared (NIR) counterparts unresolved in 
the 2MASS survey (Cutri et al. 2001). However, the res- 
olutions of the Spitzer data (2'.'3 at 8 /im, 7" at 24 /im) 
and 2MASS data (2.5") are inadequate in helping to un- 
derstand the detailed environment of these maser rings. 
Furthermore, several of the sources as seen by Spitzer 
are saturated and are so bright that they cause image 
artifacts preventing one from knowing the true nature of 
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any extended dust emission, if any. Therefore, we de- 
cided to obtain the highest spatial resolution NIR and 
MIR imaging available to explore the true nature of the 
maser rings. Using the Gemini 8-m telescopes, we ob- 
tained data with NIRI/Altair, an adaptive optics NIR 
instrument that can achieve resolutions of <150 mas at 
2 /im. MIR observations were made with T-ReCS, em- 
ploying a method which fully characterized the system 
point spread function accurately enough to allow the 
imaging data to be reliably deconvolved, achieving spa- 
tial resolutions of —150 mas at 8 /im and —250 mas at 
18 /im. 

To better understand the nature of the methanol rings, 
the main aim of our observations was to resolve the cir- 
cumstellar dust emission at NIR and MIR wavelengths 
to verify whether or not their infrared morphologies are 
consistent with the hypothesis that methanol maser rings 
are tracing circumstellar disks. Additionally, through the 
use of spectral energy distribution (SED) model fitting, 
we will learn more about the possible physical proper- 
ties and source geometries of the massive young stellar 
objects (MYSOs) associated with the maser rings. 

2. OBSERVATIONS 

2.1. Sample selection 

The targets were selected from the sample presented 
in Bartkiewicz ct al. (2009). Nine out of 31 methanol 
masers spot distributions in that study can be fit by el- 
lipses, but of those, only seven targets have four or more 
maser clumps and thus are the most convincing maser 
rings in the survey. Inspection of archived Spitzer data 
revealed that there are MIR counterparts to all seven 
maser rings. However, in four out of these seven targets 
the MIR emission in the Spitzer maps at 8 and 24 /im is 
saturated, blowing out the centers of the sources and cre- 
ating ring-like ghost structures. This makes it impossible 
to use the Spitzer data to derive any information about 
the spatial extent and/or shape of the MIR emission, 
however this did indicate that the target would certainly 
be bright enough to be observed with reasonable expo- 
sure times on less sensitive ground-based telescopes. The 
four maser ring target fields are given in Table 1. 

2.2. MIR observations 

We used the Thermal-Region Camera and Spectro- 
graph (T-ReCS) at Gemini South in queue mode (Pro- 
gram ID GS-2009B-Q-7), with data being obtained on 
several nights spanning the 2009 July 29 to 2009 August 
29 time period. Observations were taken with the Si- 
2 (A c =8.74 /im, AA=0.78 /im) and Qa (A c =18.30 /im, 
AA=1.6 /im) filters. T-ReCS employs a Raytheon 
320x240 pixel Si: As IBC array which is optimized for 
use in the 7-26 /im wavelength range. The pixel scale is 
0.089"/pixel, creating a field of view of 2 8'.' 8 x 21'.'6. Total 
exposure times were 340 s for the 8.6 /im filter and 360 s 
for the 18.3 /im filter. Standard chop-nod observing was 
used to remove sky and telescope radiative offsets. 

Flux calibration was achieved by observing the mid- 
infrared standard star HD168723 (77 Ser) at a similar air- 
mass as each of the science targets. The assumed flux 
densities were taken to be 20.67 Jy at 8.6 /im and 4.98 
Jy at 18.3 /im. These assumed flux densities were found 
by convolving the T-ReCS filter profiles with the spectral 



irradiance templates of Cohen et al. (1999). Variability 
in the measured response of the standard star from night 
to night was used to estimate the error in the flux density 
observations. From this we conclude that the measured 
8.6 /im flux densities have a 1-a error of 8% and the 18.3 
/im flux densities have an error of 15% (see Table 2). 

2.2.1. MIR Deconvolution and Resolution 

Our observations at 8.6 and 18.3 /im yielded natu- 
ral spatial resolutions of —0.29" and —0.50", respec- 
tively. Since the delivered point-spread function (PSF) 
can change due to different wind loads on the secondary 
mirror, observations of bright mid-infrared stars were 
taken immediately before and/or after each science tar- 
get observation. Furthermore the PSF can change due to 
differences in the telescope flexure as a function of tele- 
scope position, so the PSF stars were chosen to be no fur- 
ther than 1° away from the science targets. From these 
observations we believe that we can accurately character- 
ize the PSF of the science observations, and using these 
PSFs, deconvolve t he target images using the ma ximum 
likelihood method (|Richardsonlll972t iLucvl 119741 ) . That 
increases the effective resolution of our images by a factor 
of -2 (i.e., -0.15" at 8.6 /im and -0.25" at 18.3 /im). 
G23.657-00.127 and G23.389+00.185 yielded detections 
of unresolved point sources in the MIR, and thus their 
images were not deconvolved. However, G24. 634-00. 324 
and G25. 411+00. 105 displayed extended MIR emission 
and therefore we deconvolved both the 8.6 and 18.3 /an 
images of these sources. These deconvolved images com- 
pare favorably to simple unsharp masking of the origi- 
nal images, and hence all of the substructures revealed 
in the deconvolved images are believed to be real with 
high confidence. In the case of G25. 411+00. 105, the de- 
convolution was also performed on the 2.12 /im image 
to enhance the view of the complex structure observed 
there. 

2.3. NIR observations 

Three of the four maser rings (G23.389+00.185, 
G23.657-00.127, and G25.411+00.105) had bright 
enough 2MASS K band counterparts (mx < 12.9) to 
be observed with Gemini in reasonable exposure times. 
In order to image the sources in the highest spatial reso- 
lution possible at these wavelengths we employed the use 
of near-infrared adaptive-optics (AO) imaging using the 
ALTtitude conjugate Adaptive optics for the InfraRed 
(ALTAIR) in conjunction with the Near InfraRed Im- 
ager and Spectrometer (NIRI) at Gemini North. Obser- 
vations were carried out from 8 August to 16 October 
2009 in queue mode (Program ID GN-2009B-Q-44) us- 
ing the K' (A c =2.12 /im, AA=0.35 /im) filter. NIRI has 
a 1024x1024 ALADDIN InSb array and has three se- 
lectable pixel scales and fields of view. We selected the 
f/32 mode which yields a pixel scale of 0.0219"/pixel and 
a field of view of 22.4" x 22.4". Because all three sources 
lie in generally extinguished regions, there were no stars 
bright enough to guide with nearby, and the laser guide 
star (LGS) system was used. The LGS system does not 
allow the AO system to make corrections as well as it 
could with a natural guide star. Furthermore, the nights 
of our observations had less than optimal seeing. Though 
the NIRI/Altair system is capable of a top resolution of 
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0.07" at 2 /urn, our observations have resolutions in the 
0.11"-0.16" range. 

The K' band AO imaging was performed using a 3 by 
3 dither pattern with 4" offsets. The median of these 
nine images was used to produce a sky frame that was 
subtracted from each of the frames before construct- 
ing the final co-added mosaic image. The effective fi- 
nal exposure times for the three sources is 756 s for 
G23.389+00.185, 900s for G23.657-00.127, and 1080 s for 
G25. 411+00. 105. Flux calibration was achieved by ob- 
serving a nearby standard star either just before or just 
after each science observation. The standard star GSPC 
S875-C (m K > = 10.760) was used for G23.657-00.127 and 
G25.411+00.105, and FS 148 (m K > = 9.441) was used 
for G23.389+00.185. The measured flux densities of the 
calibrators varied by only 1% from night to night, and 
the derived K' band flux densities for several sources on 
each field were compared to their K values as given in 
the 2MASS Point Source Catalog. From this we estimate 
the absolute calibration errors in the measured K' band 
flux densities of our science targets to be better than 5% 
(See Table 2). 

2.4. NIR and MIR Astrometric Accuracies 

Gemini observatory has developed a technique that al- 
lows their optical telescopes to be pointed accurately at 
objects with no optical emission. The technique involves 
the use of two optical astrometric standard stars, one 
5-7' away from and the other <50" away from the non- 
optical science target. The location of the non-optical 
science target is in essence "triangulated" using these 
stars. The astrometric standard stars were chosen from 
the Second U.S. Naval Observatory CCD Astrograph 
Catalog (UCAC2) catalog, which has positional errors 
of about 20 mas for stars with m R ~ 10 — 14. Repeated 
trials using the set ups with optically bright science tar- 
gets showed that the la deviation was 0.06". 

Since it was of the highest importance that we know 
precisely where the maser rings are with respect to the 
dust emission in the MIR, we employed this astrometric 
technique to our observations. This was only done for the 
science observations at 8.6 ^m. The 18 ^m astrometry 
was achieved by registering the source positions at 18 fxm 
with the 8.6 fim emission. If the emission peak at 8.6 /im 
is actually coincident with the 18 fim peak, we estimate 
our absolute positional uncertainty is slightly worse at 
18.3 /xm (-0.10") than 8.6 /im. 

For the three sources that were also observed in the 
NIR, we employed a different astrometric technique. 
Since the NIRI fields containing the science targets also 
contained multiple NIR stars from the 2MASS Point 
Source Catalog, we were able to use their measured po- 
sitions to accurately define absolute astrometry of the 
NIRI images. A % -minimization technique was used to 
register the 2MASS fields to the NIRI fields. Then the 
offsets in RA and Dec between the final positions of the 
centroids of the NIRI sources with respect to the 2MASS 
sources were measured. The standard deviation in these 
right ascension and declination offsets yields the la error 
in the absolute astrometry of the NIR images (see Table 
I)- 

Since we are comparing the offsets of infrared emission 
from the methanol masers, we point out that the error 
in the absolute coordinates for the methanol masers is 



on the order of a few mas owing to the phase-referencing 
technique that was employed (Bartkiewicz et al. 2009). 
This means that when comparing the maser and infrared 
positions, the infrared astrometric errors quoted here for 
each source dominate in all cases. 

3. DISCUSSION AND RESULTS 

In order to fully understand how our NIR and MIR 
imaging can test the hypothesis of whether or not 
methanol maser rings trace circumstellar disks around 
massive young stellar objects, we must first state what 
we expect the IR emission from such a source to look like 
under the assumption that it is forming via disk/envelope 
accretion. We will discuss here how the observed prop- 
erties of embedded young massive stars with disks differ 
from more well-known IR observations of disks around 
low mass stars. We will then discuss the expected IR 
morphologies of a MYSO, and how that morphology dif- 
fers with angle-dependent disk geometries. We will then 
discuss the results of our high resolution IR observations 
and how they compare to the expected IR morphologies 
of the four MYSOs whose disk geometries were derived 
from the properties of their methanol maser rings. Fi- 
nally, we will use SED models to fit the multi- wavelength 
photometry of each source to derive the likely physical 
and geometrical properties of each MYSO and compare 
those results to the expected properties based on the 
masers and to the observed properties of the IR mor- 
phologies seen in the images. 

3.1. The manifestation of IR emission from a 
disk- accreting massive young stellar object 

The first infrared observations to resolve disks around 
stars were of debris disks (e.g., Telesco et al. 1988, 
Jayawardhana et al. 1998). Debris disks are circum- 
stellar disks that occur around stars with more evolved, 
(post)planet building disks. Since one can directly view 
a debris disk in its infrared dust emission it appears as 
an elongated ellipsoidal structures whose "flatness" is 
a function of the viewing angle of the disk. However, 
at earlier stages of star formation, such circumstellar 
disks are likely to be actively accreting on to their par- 
ent stars, and are thus much more massive and dense. 
In fact, accretion disks around very young stars began 
to be found, not by their direct emission, but by their 
scattered and reflected emission off of their upper and 
lower disk surfaces. Originally discovered in the opti- 
cal (e.g., McCaughrean & O'Dell 1996), these "silhou- 
ette" disks demonstrated that if a disk is large and dense 
enough, it could be optically thick even at infrared wave- 
lengths (e.g., McCaughrean et al. 1998, Cotera et al. 
2001). These sources are typified by a "dark lane" de- 
marcating the disk itself, between two "lobes" of scat- 
tered/reradiated emission demarcating the upper and 
lower flared disk surfaces. 

However, all of these observations were of disks around 
low-mass stars. The first claims of infrared detections of 
circumstellar disks around young massive stars (Stcck- 
lum & Kaufl 1998, De Buizer et al. 2000) were made, in 
part, on the assumption that the disk would appear as 
an elongated structure in its infrared dust emission, sim- 
ilar to what was seen with debris disks around low-mass 
stars. However, higher resolution follow-up observations 
(De Buizer et al. 2002) proved that these were not disks, 
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and it was soon after realized that disks accreting onto 
young massive (proto)stars would be optically thick in 
the infrared as well, and would be even more dense and 
massive than those accreting onto low-mass protostars. 
Furthermore, in order to have a large enough reservoir 
of material to accrete to such high masses, massive stars 
must form in the densest parts of giant molecular clouds, 
which are extremely obscured environments. Moreover, 
the earliest stages of massive star formation arc deeply 
self-embedded; the (proto)stars are surrounded by mas- 
sive and dense accretion envelopes, as well as the afore- 
mentioned circumstellar disks, which are believed to be 
large, thick, and flared. 

Due to the high amount of obscuration, one would have 
no hope of observing these stages of massive star forma- 
tion at even MIR wavelengths if not for their bipolar out- 
flows. Disk accretion is accompanied by outflow, and it is 
this outflow that punctures holes through the obscuring 
material surrounding a forming massive star. The out- 
flow axis is oriented perpendicular to the disk, and the 
outflow starts out narrow and collimatcd. At such early 
stages of formation, we may only detect the presence of 
the massive young stellar object at infrared wavelengths 
if we are lucky enough to have a line of sight looking 
down the outflow axis. We could then see into the enve- 
lope, and view the scattered/reprocessed dust emission 
off of the cavity walls, and if the angle is just right, we 
may see down into the central disk or (proto)star. Such 
a chance alignment is rare, and it is in part due to this 
geometrical effect that detecting massive young stars in 



sonable heating argument 2 . Furthermore the presence of 
methanol masers is thought to be a signpost of the earli- 
est stages of MYSO formation (e.g., Breen et al. 2010). 
Additionally, the youngest MYSOs are known to display 
no cm radio continuum, and none of our sources were 
found to be associated with 8.4 GHz radio continuum 
emission (Bartkiewicz et al. 2009). This means that if 
the methanol masers in our sample are within circumstel- 
lar disks in these sources, then they likely trace a rather 
young stage of the massive star formation process where 
the disk is an active accretion disk, likely flared, sur- 
rounded by a thick accretion envelope. Figure 1 shows 
a simplistic toy model of this, where the outflow cavi- 
ties are the source of the IR emission from the MYSO. 
The disk and the envelope are assumed to have a uni- 
form opacity (i.e. no formal radiative transfer calcula- 
tions are implemented), and the cavities are assumed to 
have uniform emission. Of course these conditions are 
unrealistic, but show to first order the bulk observable 
properties that are the result of the more formal mod- 
eling of Alvarez et al. (2004) and Zhang & Tan (2011), 
and thus are adequate for demonstration purposes. 

The toy modeling in Figure 1 shows that there are ex- 
pected spatial and morphological relationships that can 
be tested between what one sees in infrared emission 
(blue and red areas in Figure 1) when looking toward 
methanol masers rings (green circles) under the assump- 
tion that the rings are in disks and that they are asso- 
ciated with an early and embedded accretion phase of 
massive star formation: 



the infrared is so difficult. However, over time the young 
stellar object evolves and the outflow angle widens (SlnScenario 1: The more circular the maser ring, the more face- 
& Adams 1987, Beuther & Shepherd 2005). With the on the disk, meaning one would be looking right 

widened outflow, comes a large range of angles that al- down the blue-shifted outflow axis into the outflow 

low for a higher probability of detection in the infrared. cavity. Therefore, the infrared emission should be 

At some point the angle is so wide, that the distinction unresolved or circularly symmetric with the peak 

between what is the outflow cavity surface and what is of infrared emission coincident with the maser ring 

the surface of the flared disk becomes a matter of seman- center; 
tics. 

This notional sequence of events has been supported cenario 2: The more highly elliptical the maser ring, the closer 



by several studies in the recent decade. For instance, ob- 
servations of the earliest stages of massive star formation 
with collimated outflow cavities emitting brightly in their 
mid-infrared continuum emission were first identified by 
Dc Buizer (2006). The first claim of a candidate infrared 
"silhouette disk" around a massive star was made by 
Chini et al. (2004) demonstrating the later stages where 
the outflow has widened considerably (though whether 
the mass of the central object is high mass is the subject 
of debate, e.g. Sako et al. 2005). In addition to these ob- 



to edge-on will be the disk. In the edge-on case the 
orientation makes it difficult to detect the MYSO 
in the infrared, unless the envelope is not opti- 
cally thick to the IR wavelength in question. If de- 
tectable in the infrared, we would expect for mod- 
erately to highly elliptical maser rings to see some- 
thing more like a silhouette disk in the infrared, 
where the maser ring would lie between two in- 
frared bright sources (the outflow cavities), in the 
"dark lane" of the optically thick disk; 



servations, MYSO radiative transfer models of Alvarez et 

al. (2004) and Zhang & Tan (2011) also reproduce thest? cenario 3: For an intermediate ellipticity maser ring, one 



geometries; the first showing the results of the earliest 
stages of collimated outflow cavities, and the second for 
more open-angled outflows. 

3.2. Testing the disk hypothesis of maser rings with IR 

imaging 

Because methanol masers are believed to be pumped 
by mid-infrared photons from dust (Cragg, Sobolev, & 
Godfrey 2002), if the masers are arising from gas in a 
dusty circumstellar disk, they would have to trace the 
inner few hundred AU radius at most, given any rea- 



would expect an intermediate disk inclination. 
In this case we would expect to either only see 
the blue-shifted outflow cavity, or bright emission 
from the blue-shifted cavity with a fainter com- 
ponent seen from the red-shifted cavity (with the 
brightness of the red-shifted cavity dependent on 
disk/cnvclopc extinction, outflow opening angle, 
and disk inclination). In this case, the maser ring 

2 Remember however, the flaring of the disk and the dense ac- 
cretion envelope may prevent one from viewing the mid-infrared 
emission coming directly from the disk for most, if not all, viewing 
angles. 
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center will be slightly offset from the infrared emis- 
sion center in a direction given by the outflow axis. 

3.3. Results from the high resolution IR imaging 

Figure 2 sho ws the four 6.7 GHz meth anol maser rings 
in our sample (Bartkic wicz et aT1l2009f) . Next to each 
plot of the maser rings the implied disk properties are 
listed based on fits to the maser rings. Also shown are 
the toy models created from those same properties, which 
demonstrate the gross expected infrared emission mor- 
phology for each source. In Figures 3-6 we present maps 
of NIR emission towards three targets (G23. 389+00. 185, 
G23.657-00.127 and G25.411+00.105) and MIR emis- 
sion at 8.6 and 18.3 fim towards all four targets. In three 
out of the four cases we found that the MIR emission at 
8.6 /xm is either extended or from multiple components. 
Furthermore, in all three cases where we also observed 
the maser rings in the NIR at 2.12 /jm, we were able 
to resolve the emission with the high spatial resolution 
afforded by adaptive optics. In Table 1 we list the co- 
ordinates of component peaks detected in the NIR and 
MIR 3 , and in Table 2 we list the integrated flux densi- 
ties of the components in the field of view at all available 
bands. 

Below we describe results of the infrared observations 
of each of the four methanol maser rings relevant to the 
aims of this project. The emphasis here is on the whether 
or not the NIR and MIR emissions seen are consistent 
with the scenarios developed in Section 3.2 and toy mod- 
els shown in Figure 2. There are a few other interesting 
details that emerged from the data that are not impor- 
tant to the discussion here; these are summarized in the 
Appendix. 

3.3.1. G23. 389+00. 185 

This source has an inferred disk inclination of 54°, so 
in this case we would expect to see infrared emission 
distributed as discussed in Scenario 3. 

We find a bright unresolved MIR source (source 1) with 
a peak close (272 mas) to the maser position, and a sig- 
nificantly weaker source (source 2) 2" (^9000 AU) to 
the northeast of this main source in both bands (Figure 
3). This second source is far enough away from source 
1 that they are likely not two IR-bright cavities from a 
bipolar outflow of a single MYSO, and are more likely to 
be separate (proto)stellar sources. Further evidence of 
this comes from the fact that source 2 appears to have 
its own jet emanating from it, as seen in the NIR (see 
Figure 4; details are given in the appendix). 

The NIR observations also reveal both sources, with 
source 1 being a resolved, circularly symmetric source 
(although this particular NIR data had elongated image 
quality; see Figure 4). The peaks of source 1 and 2 in 
the NIR coincide with the peaks in the 8.6 /jm image to 
better than 60 mas in both cases, and thus the MIR/NIR 
peaks are believed to be co-spatial. This demonstrates 
the effectiveness of both of the two independent astro- 

3 The naming convention of the IR sources is from the IAU 
recommendation, which take the form Gll.lll±bb.bbb:BDS12 #, 
where the # represents the component number. These names are 
listed in the tables for formality, but in the figures and when dis- 
cussing components within a particular field are often abbreviated 
to just the component number alone. 



metric techniques, i.e., the MIR astrometry technique de- 
scribed in the observations section and the ability to tri- 
angulate NIR source positions from other 2MASS sources 
in the field. It also leads to the conclusion that the offset 
between the maser ring center and infrared peak is sig- 
nificant in this case. In fact, because there was very good 
spatial agreement between the positions of the 2MASS 
sources and the NIRI sources on the field in declination, 
the astrometric technique (described in Section 2.4) led 
to very small errors along that direction for this source 
(15 mas). Since the offset of the maser ring center and 
the NIR emission center is mostly in declination, that 
offset is 13(7. 

Moreover, the angle on the sky that the IR peak of 
source 1 makes with respect to the maser ring center (p. a. 
= +10°) is quite different than the disk/outflow axis as 
inferred from the maser ring geometry (p. a. = —45°), 
inconsistent with Scenario 3 and the toy model shown 
in Figure 2 for this source. However, in Scenario 2 the 
maser center and IR peak don't have to be coincident, 
but simply lie along the disk/outflow axis. The closest 
point along this axis to the IR peak is indeed closer (210 
mas), but would still be significantly more than allowed 
by the combined NIR/MIR astrometry errors (>4.5<r). 

The only outflow imaging observations conducted to- 
ward this region were the low resolution (~27") obser- 
vations of HCO+(3-2) by Schenck et al. (2011). Both 
IR sources 1 and 2 lie on the southeastern edge of the 
extended region of integrated HCO+ emission, with the 
peak located ~20" to the northeast of source 2. Given 
the coarse resolution it is not clear if the integrated 
HCO+ emission is associated with IR source 1 or 2, 
and with no channel maps it is unclear if there are well- 
defined red and blue shifted lobes. It is plausible given 
the location of the the integrated HCO+ peak that it 
traces the outflow emission associated with the counter- 
jet to the NIR jet seen coming from source 2. If this were 
the case, the HCO+ would be unrelated to the masers 
associated with IR source 1. If the HCO+ emission is 
tracing an outflow lobe from the YSO associated with the 
masers, the angle the HCO+ peak makes with respect to 
the maser ring center would imply a outflow direction 
perpendicular to what is expected from the maser ring 
geometry. 

In summary, the infrared (and HCO+) emission from 
this source does not seem to be consistent with the as- 
sumption that the maser ring traces a circumstellar disk. 

3.3.2. G23.657~00.127 

For this source, the maser ring is almost circular with 
an infe rred disk inclination of o nly 16° from disk model 
fitting (jBartkiewicz et aT1l2005l ). It is therefore expected 
that the IR emission will manifest itself as described in 
Scenario 1. 

Consistent with that scenario, emission from a bright, 
unresolved source (source 1) is seen in the MIR maps 
of this maser field (Figure 4). We find that the peak of 
8.6 ^m emission and the peak of the 2.12 /^m emission are 
coincident with one another to within 0.11", or 0.8(1.8)<r 
of the NIR(MIR) astrometry error, and thus are likely 
to be co-spatial. Consistent with Scenario 1 and the 
toy model for this source, the center of methanol ring 
is nearly coincident with the NIR(MIR) peak to within 
1.0(2.5)ct. 



6 



Dc Buizer, Bartkiewicz, & Szymczak 



However, contrary to what is expected in Scenario 1, 
at 2.12 fj,m source 1 is resolved into a extended triangu- 
lar fan of emission. The top panel of Figure 2 demon- 
strates that this morphology is real, as there are two 
other moderately bright sources on the NIRI field that 
are point-like. While ideally an outflow cavity viewed 
close to pole-on should have circular symmetry, a non- 
circular emission morphology is not in itself surprising; 
in reality, environmental conditions will likely decrease 
this level of circularity. However for this source, the NIR 
emission is highly flattened along two sides (see Figure 
4), strongly implying the presence of a rather straight 
barriers of some sort. At present no outflow or molecu- 
lar imaging of sufficient angular resolution of this object 
exist to help with the interpretation. However, since the 
extended fan of emission is seen only in the NIR, the sim- 
plest explanation is that it is likely scattered and/or re- 
flected emission off the walls of the outflow cavity. There- 
fore, given the NIR morphology, this source seems more 
consistent with a much more inclined disk/outflow than 
what one would infer from the maser ring geometry. 

Therefore, the observed infrared morphologies as a 
whole are not fully consistent with the hypothesis that 
the maser ring arises from a face-on circumstcllar disk 
around this source. 

3.3.3. G24.634- 00.324 

In the case of this maser ring, whose inferred inclina- 
tion from a fitted disk model is 71°, we expect to see an 
infrared morphology close to that described in Scenario 
2. 

We do indeed detect a bi-lobed structure in the MIR 
(Figure 5), with two bright peaks at 18.3 and 8.6 ^m 
separated by ~1" (no NIR data were obtained for this 
source). This could be interpreted as a silhouette disk- 
like morphology, where the dark lane between the sources 
would be the disk midplane. Indeed, as can be seen 
in the natural resolution 8.6 /im image in Figure 5, the 
methanol maser ring appears to be co-located with the 
dark lane to within 1.2er of the astrometric uncertainty. 
The size of the dark lane extends ^1500 AU in diameter, 
which would be large for a inner accretion disk, but is 
comparable to what one might expect for a larger-scale 
outer flared disk or torus (Bcuther et al. 2009). The de- 
convolved 8.6 (im image shows the sides of sources 1 and 
2 along this dark lane to be very flat, which may be ex- 
pected if a disk is present, but which would be unusual if 
these two sources were individual MYSOs. Lending fur- 
ther credence to the silhouette disk-like interpretation is 
the wavelength dependent nature of the MIR emission. 
Both sources 1 and 2 grow brighter with wavelength (i.e. 
18.3 /im compared to 8.6 /xm), but source 2 increases in 
brightest more and becomes more extended. This would 
be expected if source 1 is the side of the disk/outflow 
facing us and source 2 is the back side partially obscured 
by the extended and flared disk. The disk becomes less 
opaque at longer wavelengths, thus would influence the 
brightening of source 2 more than source 1 as it is ob- 
served at longer and longer wavelengths (see Zhang & 
Tan 2011). The 1.2 mm continuum maps of this source 
(Rygl et al. 2010) show a that the cool dust emission in 
this area peaks ~6" to the northeast of the maser ring po- 
sition. However, the mm continuum source is elongated 
at a position angle of ~60°, which is approximately the 



position angle of the MIR "dark lane" , and extends over 
~30". Optimistically, one could claim that this is con- 
sistent with the "silhouette disk" interpretation of the 
IR emission, however because the resolution of the mm 
maps is only 10.5", it is not known if the elongation in the 
mm emission is due to a large torus or simply unresolved 
multiple dust condensations. 

However, if this is indeed a silhouette disk-like mor- 
phology, the angle of the outflow/disk axis given by the 
MIR morphology is 90° from the axis inferred from the 
methanol maser ring from the best fitted ellipse. There- 
fore, the IR morphology and mm dust emission is not 
consistent with the toy model for this source (as given in 
Figure 5) and therefore inconsistent with the maser ring 
tracing a circumstellar disk. 

Though we feel that the "silhouette disk" disk inter- 
pretation is the best explanation for the IR emission 
given the limited data at hand, it is however plausible 
that the two sources of IR emission are separate, indi- 
vidual MYSOs. If this were the case then the maser ring 
lies closer to, and would be more likely associated with, 
source 2. Some credence to this idea comes from the 
presence of source 3. This source disappears at longer 
wavelengths, signifying that the source is hotter and/or 
less embedded star. This could lead to the idea that all 
three sources are a tight cluster of YSOs with source 2 
being the closest YSO to the masers. However, the maser 
ring center lies more than 3.2<r (of the astrometric un- 
certainty of the MIR astrometry) away, i.e. much farther 
away than the "dark lane" . Ignoring this, if source 2 is an 
individual MYSO associated with the masers, its emis- 
sion could be interpreted as the face of a disk/outflow 
cavity facing towards us with the opposite face/cavity 
completely obscured from view, in other words the mor- 
phology described in Scenario 3. Looking at Figure 5, 
the MIR peak of source 2 does not appear to lie along 
the disk/ouflow axis as inferred from the maser ring ge- 
ometry. However, the maser ring center is only 2Aa away 
from the closest point along this axis. Thus this scenario 
could be plausible (however there is even further evidence 
against this; see Section 3.4), and could be consistent 
with the maser ring residing in a disk, if the sources of 
emission seen are individual MYSOs. 

Source 3 may also be a foreground or background 
source, but the 2MASS and Spitzer IRAC images do not 
have the resolution to confirm this. 

On a final note, since this source has only four maser 
groups fit by an ellipse (see Figure 2), its categorization 
as a maser ring is the most tenuous of the four sources 
studied. 

3.3.4. G25. 411+00. 105 

This maser ring has an inferred disk inclination of 47° 
from disk modeling, and therefore is expected to have an 
IR emission morphology as described in Scenario 3. 

The 2.12 jum map reveals a complex structure of NIR 
emission. Using a bright point source on the field, the 
2.12 fim image was deconvolved to further enhance the 
resolution (to ~80 mas). This is shown in Figure 6. The 
NIR peak labeled 3 corresponds most closely to the loca- 
tion of the maser ring center, with only a 1.9a (200 mas) 
offset, but this offset is in the direction expected from 
the toy model for this source (Figure 6). Consistent with 
Scenario 2 and the toy model, the bulk of the NIR emis- 
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sion seems to be distributed along the disk/outflow axis 
as defined by the fit to the maser ring (i.e., p. a. = 180°). 

However, in the MIR images we see quite a different 
morphology (Figure 6). At both 8.6 and 18.3 fi,m, there is 
a well-defined peak, with an extended tongue of emission 
to the northwest. This main MIR peak is only 70 mas 
away from the NIR peak labeled source 3 and therefore 
assumed to be co-spatial. 

The tongue of emission seen to the northwest of the 
peak in the MIR, is coincident with the NIR sources 6 
and 7. In fact there are local weak peaks in the de- 
convolved MIR images that appear to be co-spatial with 
these NIR peaks. This MIR emission is not what would 
be expected for the morphology described in Scenario 3 
or the toy model for this source, and therefore it does 
not seem to support the idea that the maser ring is in 
a circumstellar disk around source 3. However, further 
observations may be needed to determine if sources 6 
and 7 are separate young stellar sources, or in any other 
way emission unrelated to the young stellar object with 
the maser ring. If they were found to be unrelated, this 
would give credence to the hypothesis that the masers 
here may be arising within a circumstellar disk. 

Outflow observations do exist for this source but un- 
fortunately not at sufficient resolution to clarify the sit- 
uation. The outflow is seen best in the CO (2-1) maps 
of Beuther et al. (2002), which shows blue-shifted and 
red-shifted lobes with peaks situated at a position an- 
gle of ^145°. The midpoint between these two peaks is 
~15" to the northeast of the maser location, and there- 
fore any YSO directly associated with the masers cannot 
be the main source driving the outflow. However, the 
CO contours have secondary peaks in both outflow lobes 
situated about ~15" to the west, with the blue-shifted 
CO outflow secondary peak coincident with the maser 
location and the MIR emission. There may therefore be 
two outflows located near each other in this region, but 
unresolved in the CO maps, one of which may be asso- 
ciated with the masers. The integrated HCO+ emission 
maps by Schenck et al. (2011), appear to only trace 
the blue-shifted lobe of the outflow, but interpretation 
is difficult with the limited resolution (~27") of those 
maps. Further complicating the outflow picture is the 
fact that there is another bright NIR/MIR source 6" to 
the north of the maser ring, which is likely a YSO unre- 
lated to the masers. The position of this source is close 
to the midpoint between the secondary red-shifted CO 
peak and secondary blue-shifted peak, and could poten- 
tially be the source of the outflow, rather than the YSO 
associated with the maser ring. 

In conclusion, while some evidence appears to support 
the disk hypothesis for this source, the MIR emission is 
not at all consistent with that picture. 

3.4. Results from SED model fitting 

The qualitative comparisons made so far between the 
expected IR morphologies and the observed IR mor- 
phologies are a good first-order test of the maser ring 
hypothesis. The next step is to see if the IR emission ob- 
served for each source quantitatively makes sense with re- 
spect to the maser ring hypothesis. A first cut at this can 
be done with existing YSO spectral energy distribution 
(SED) models, insofar as the change of the integrated 
flux density of an object as a function of wavelength will 



vary substantially as a function of MYSO geometries, in- 
cluding disk inclination and outflow cavity opening angle. 
Through the use of our observed flux densities and those 
available from data archives, we should be able to con- 
strain the modeling enough to at least rule out certain 
geometries. 

In order to derive and interpret the SEDs, we ap- 
plied the mode l of yo ung stellar objects presented by 
iRobitaille et al.l (|2007|) 4 . First we searched for coun- 
terparts of the four targets in the following catalogues: 
2MASS, Spitzer IRAC, MSX, and IRAS. For all four 
maser sources, the detected IRAS fluxes were used as up- 
per limits in the SED models, due to the extremely coarse 
resolution of the IRAS data (~1') and thus their poten- 
tial for contamination. The 2MASS, IRAC, and MSX 
data were also taken with coarser resolutions than our 
Gemini NIR and MIR images (as mentioned in Sect. 1), 
however those same Gemini data were used to determine 
if these archival data would be used as upper limits or 
not. The source IR fluxes and distances that we used in 
the SED modeling are listed in Table 3, and data used 
as upper limits are marked. 

The model by IRobitaille et al.l (|2007t ) is based on the 
computed 20,071 YSO radiation transfer models at ten 
different viewing angles that span a large range of evo- 
lutionary stages and stellar masses. All of the main fea- 
tures of our toy model are formally integrated into these 
models; they use pre-main-sequence stars with different 
combinations of axis-symmetric circumstellar disks, in- 
falling flattened envelopes, and outflow cavities under 
the assumption that stars form via accretion through the 
disk and envelope. The fitting procedure interpolates the 
model fluxes to the apertures used in measurements (the 
source sizes) scaling them to a given distance range of a 
source. It also allows the interstellar extinction, Ay, to 
be a free parameter in the fitting process, which we chose 
to be 0<Ay<100mag. 

We present in Figure 7 the model fits of the SEDs. 
In Table 4 we list the derived physical properties of the 
modeled MYSOs, including stellar mass (-M*), disk in- 
clination (idisk), disk outer radius (Rdisk), outflow cav- 
ity opening angle (6 'cavity) i an d total luminosity (L to t)- 
Several other properties are also given by the SED fits, 
however many of them are degenerate or not as impor- 
tant for the comparisons made here. For each source we 
examine the best 10 SED fits based on their \ 2 values 
and list the best fitted values, the modes and medians for 
each parameter, and also give the minimum and maxi- 
mum values of each parameter. It can be seen in Table 
4 that the central young stellar objects are indeed mas- 
sive (8-15 Mq for the best fits), as expected for sources 
exciting methanol maser emission. 

In Table 4, we also present the disk inclination and 
radius as derived from the maser ring geometry for each 
source. These can be directly compared to the MYSO 
geometries derived from the SED modeling to check for 
consistencies/inconsistencies of the maser disk hypothe- 
sis or our interpretations made so far based on the qual- 
itative morphologies of the IR emission. 

In the case of G23. 389+00. 185, we note a good agree- 
ment between the photometry from the high resolution 
MIR and NIR images and with the archive data. This 

4 The tool available via http:/ /caravan. astro. wisc.edu/protostars, 
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translates into a fairly well constrained fit to the fluxes 
(Figure 7). This also can be seen in the similarities be- 
tween the best fit, the mode fit, and the median fit in 
Table 4. The problem with this fit is that it seems to be 
quite different from what would be expected if the maser 
ring were tracing a disk. The maser ring geometry im- 
plies a disk inclination of 54° and a radius of 425 AU, 
however the models show that the most likely fits come 
from a more face-on geometry with a disk inclination of 
18° and a radius less than 100 AU (and in most cases no 
disk at all). This adds further evidence to our claim from 
the discussion of the IR morphologies that the maser ring 
is not consistent with the disk hypothesis. 

For G23. 657— 00.127, we again note a good agreement 
between the Gemini MIR and NIR photometry with 
archive data, which again translates into relatively well- 
constrained SED results. However, the results do ex- 
hibit a bimodal nature, with two groupings of tightly 
constrained parameter values that can be separated by 
disk inclination. For this source, therefore, we show the 
results of two top ten lists in Table 4. The first group of 
results are for the more face-on disk inclination, and the 
second for the more inclined disk results. It can be seen 
from the statistics of these two groups that there is very 
little spread in properties once the results are separated 
by the two disk inclinations. This is also evident in the 
SED plot for this source (Figure 7). 

Interestingly, this dichotomy in the SED model results 
parallels the results from Section 3.3.1 based solely on 
the IR morphologies for this source. From the maser ring 
geometry alone, one would expect a nearly face on disk 
(i = 16°). Recall that the MIR morphology by itself is 
consistent with a face-on geometry. Consistent with that 
picture, we find here that the best SED model is for a 
disk with a inclination of 18°. Conversely, however, recall 
that the NIR emission morphology is fan-like, and looks 
more like a more like a wide outflow cavity inclined to 
the line of sight. We see here that this geometry is more 
consistent with the majority of the top 20 SED modeling 
results which yield a disk with a high inclination (41- 
76°). Moreover, the SED results for an inclination of 18° 
seem to create further problems for the maser ring hy- 
pothesis for this source, because for all of these fits (save 
the best fit) require diskless systems. In other words, 
the SED models that predict a low system inclination 
also predict an extremely small or diskless system, and 
therefore seem to contradict the maser disk hypothesis 
for this source. 

For G24. 634-00. 324 the data did not constrain the 
models very well. Under the assumption that the IR 
morphology is due to an edge-on system, we used the 
combined flux of all emission on the Gemini fields in the 
SED modeling. The results all fell into three categories, 
the first being disk-only systems, the second being nearly 
face-on geometries, and the third being the rest of the 
results. We can, from the morphology in the high reso- 
lution IR images alone, ignore all results from fits that 
imply a face-on geometry. Second, as we have discussed 
in Section 3.2, we have good reason to believe that these 
sources are very young MYSOs still embedded in their 
accretion envelopes, so we also ignore the SED results 
that are consistent with disk-only systems. The top 10 
remaining fits are shown in Table 4 and Figure 7. The 
parameters are consistent with a relatively inclined sys- 



tem (i = 41 — 76°) and therefore lend some credence to 
the belief that the MIR emission seen in the Gemini im- 
ages comes from a "silhouette-like" disk. This however, 
would also confirm that the maser ring is elongated at 
the wrong angle to be tracing the disk in the "dark lane" 
as seen in the IR. If we just use the Gemini fluxes for 
source 2 and use the Spitzer and MSX data as upper lim- 
its (since they do not resolve source 1 and 2) the SED fits 
imply a central stellar mass less than 8 M©. This seems 
to again strengthen the case for the IR emission mor- 
phology being from a MIR "silhouette-like" disk (rather 
than a cluster of three individual YSOs) and that the IR 
images are inconsistent with the maser ring hypothesis. 

In the case of G25. 411+00. 105, we assume that all of 
the flux in the Gemini images is from the outflow cavity 
of a single source. However, because of the close prox- 
imity of another bright IR source ~6" to the north, the 
~18" resolution MSX data have been used as upper lim- 
its (Table 3) for these fits. The results from the fits 
are mixed; some parameters seem well-constrained (i.e. 
idisk and M*), but others are not. Interestingly, all fits 
yield an inclined disk geometry that is consistent with 
the maser ring inclination. However, the derived cavity 
opening angle (9 cavity = 6°; see Table 4) may be too 
narrow to explain the IR emission as all coming from the 
blue-shifted outflow cavity. Furthermore, the maser ring 
is 550 AU in radius, and the median disk size of the fits is 
only a fourth of that, with the best fit being only 12 AU. 
This casts further doubt on the interpretation that the 
maser ring is in a circumstellar disk around this source. 

3.5. Do methanol maser rings trace circumstellar disks? 

Of the four methanol maser rings we observed in the 
infrared, none convincingly displayed the characteristics 
expected if the maser rings are indeed tracing circumstel- 
lar disks. In the case of G23. 657-00. 127, the maser ring 
is nearly circular, indicating that we should be looking 
almost straight down the outflow cavity (i.e., Scenario 
1). However, we see that the 2 /im morphology is fan- 
shaped as one would expect for the reflected light from a 
outflow cavity with a much more inclined geometry. The 
majority of SED model fits to the infrared photometry 
of this source seem to corroborate this more inclined disk 
geometry. In the case of G24. 634-00. 324, the maser ring 
is highly elliptical, indicating that we should be seeing a 
more edge-on disk geometry, and perhaps a silhouette- 
like infrared morphology where the masers reside in the 
"dark lane" of the infrared emission (i.e., Scenario 3). 
Consistent with the results from the SED modeling, we 
do indeed see this type of morphology in the MIR how- 
ever the angle of the "silhouette disk" is almost 90° from 
the expected disk/outflow axis as given by the maser ring 
geometry. For G23. 389+00. 185, the maser ring is slightly 
elliptical, indicating that we will likely only see the disk 
face/outflow cavity facing towards us (i.e. Scenario 2). 
However, the maser ring center is significantly offset from 
the IR peak and at an orientation that is not consis- 
tent with emission from a disk/outflow as derived from 
the maser ring geometry. Finally, for G25. 411+00. 105 
which also has a maser ring with a slight ellipticity (i.e. 
Scenario 2), has NIR emission consistent with scatter- 
ing/reflection off of a extended outflow cavity and at the 
right orientation as inferred from the maser ring geome- 
try. However, the MIR emission is elongated at an angle 
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130°diffcrcnt from this, and has a morphology that casts 
doubt on the validity of the conclusion that would be 
drawn from just the NIR emission and maser ring geom- 
etry alone. Furthermore, SED models have a hard time 
reproducing the large disk size implied by the maser ring 
geometry. 

In summary, none of the four sources observed 
give compelling evidence in support of the hypothesis 
that maser rings may arise from circumstellar disks. 
Moreover, the two clearest examples of the maser 
ring morphology in our sample, G23. 657-00. 127 and 
G23. 389+00. 185, seem to have fairly good observational 
evidence that contradicts such claims. Therefore the con- 
clusion from these observations is that the infrared emis- 
sion from these sources does not seem to support the 
scenario where methanol maser rings trace circumstellar 
disks around young massive stars. 

If maser rings are not tracing disks, what else could 
explain these ring structures? As discussed in the intro- 
duction, the scenarios that have been proposed to gener- 
ally explain 6.7 GHz methanol maser emission are disks, 
outflow and shocks. If the maser rings are not in disks, 
then they may be coming from shocks within outflows. 
Ellipses are created geometrically from the cross-section 
of a cone and a plane. Given the conical geometry of 
the outflow cavity, a quasi-planar accretion shock ema- 
nating from the (proto)star/disk could travel down the 
outflow cavity walls and cause a geometrically ring-like 
shock. This shock would liberate methanol from the 
dust, which could then be radiatively pumped to pro- 
duce the maser emission. Though in principal this is 
a possibility, our data here cannot adequately address 
this hypothesis. The outflow cavities of MYSOs can be 
heated by MIR emission out to 10s of thousands of AU 
(De Buizer 2006), and therefore the ring produced in this 
way could be located anywhere along the outflow cavity. 
Furthermore, due to possible differences in density on 
opposite sides of an outflow cavity's walls, the propaga- 
tion of shock down the cavity can skew the ring so that 
the minor axis of the ellipse does not point back to the 
central (proto)star/disk. This means that the location 
of the IR emission with respect to the maser ring center 
and with respect to the maser ring geometry are not well- 
constrained, unlike with the disk hypothesis. Therefore, 
with only the data we have here we cannot truly test the 
outflow shock hypothesis. 

It seems the obvious way forward to solve the mys- 
tery of what causes the maser rings will be to employ 
ALMA to directly image the disks at high spatial resolu- 
tion around these sources. Given directly observed disk 
geometries and their location with respect to the maser 
emission, one may finally be able to disentangle the true 
nature of methanol maser rings. 

4. CONCLUSIONS 
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The goals of this paper were to resolve the NIR and 
MIR emission at the locations of four methanol maser 
rings to verify whether or not their IR morphologies are 
consistent with the hypothesis that methanol maser rings 
are tracing circumstellar disks, and to derive physical 
properties and geometries for their associated massive 
young stellar objects through SED model fitting. 

In this article we discussed the assumed distribution 
of circumstellar material around such young and mas- 
sive accreting (proto)stars, and what infrared emission 
geometries would be expected for different disk/outflow 
orientations. For the four targets we observed, we com- 
pared the expected infrared geometries (as inferred from 
the properties of the maser rings) to actual high spatial 
resolution near-infrared and mid-infrared images. We 
find that the observed infrared emission geometries are 
not consistent with the hypothesis that the masers are 
residing in circumstellar disks. 

Using SED model fitting, we found that the emission 
from the infrared counterparts for all methanol masers 
ring distributions arc indeed consistent with massive 
young stellar objects with masses above 8M . Further- 
more, we find that in most cases the geometries allowed 
by the SED model fits corroborate the negative results 
from the observed infrared morphologies, casting further 
doubt on the hypothesis that the methanol maser rings 
in these four cases arise from within circumstellar disks. 
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APPENDIX 

ADDITIONAL INFORMATION ON SOURCES 

During the study of the sources in this article, interesting results were gathered that do not have any direct relation- 
ship to the main goals of the paper. However, this information is likely to be of use to those who study these sources 
in-depth. Here below we summarize auxiliary information we have for two sources in this study. 

G23. 389+00. 185. The MIR emission towards this source is the brightest among five targets. We note this 
methanol ring was also the brightest maser in the selected sample of methanol sources with S p =21.55 Jybeam -1 
(Bartkiewicz ct al. 2009). We also found the 22 GHz water masers towards this source (the only detection of the four 
in the sample) w ith a distribution along a position angle of 45° from the major axis of the fitted ellipse to the 6.7 GHz 
methanol spots (Bart kiewicz et al. 2011). However, their position uncertainties are of order 0'.'15. 

Interestingly, there is a bright jet of emission seen in this image, that can be better seen in a unsharp mask of the 
data (see Figure 2). It appears to point s traight back to source 2. Jets from outflowing massive young stellar objects 
have be observed in infrared continuum (Dc Buizer ct al. 2006). However, this jet-like structure is not seen in the 
MIR, thus indicating that if it is a jet it must be line emission, not continuum. The K' filter of NIRI is centered on 
the 2.122 fim v = — 1 S(l) line of H2, a well-known outflow indicator. This therefore likely a jet coming from source 
2. 

G25. 411+00. 105. The brightest NIR peak, source 1, has no MIR counterpart, and there are many other peaks 
identified in the NIR that also have no MIR counterparts. This means that, like the fan-like extended structure of 
G23. 657— 00.127, the NIR emission could be scattered emission, or alternatively it may line emission, as discussed for 
G23.389+00.185 above. 
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TABLE 1 

Maser. Fields and the Positions of Detected Infrared Sources 



NIR (2.12 /mi) MIR (8.6 ^m) 



Peak Position (J2000) la Error Peak Position (J2000) la Error 

Maser Field : IR Source a 8 (Act,A<5) a <5 (Aa,A5) 









(h m s) 


(° ' 


") 


(mas, mas) 


(h m s) 


(° ' ") 




(mas, mas) 


G23.389+00 


185:BDS12 


1 


18 33 14.323 


-08 23 


57.30 


155,15 


18 33 14.327 


-08 23 57 


29 


60,60 




BDS12 


2 


18 33 14.417 


-08 23 


55.76 


155,15 


18 33 14.417 


-08 23 55 


78 


60,60 


G23.657-00 


127:BDS12 


1 


18 34 51.561 


-08 18 


21.52 


106,157 


18 34 51.556 


-08 18 21 


51 


60,60 


G24.634-00 


324:BDS12 


1 










18 37 22.683 


-07 31 41 


58 


60,60 




BDS12 


2 










18 37 22.713 


-07 31 42 


33 


60,60 




BDS12 


3 










18 37 22.659 


-07 31 42 


68 


60,60 


G25.411+00 


105:BDS12 


1 


18 37 16.920 


-06 38 


31.00 


114,96 










BDS12 


2 


18 37 16.908 


-06 38 


30.74 


114,96 












BDS12 


3 


18 37 16.917 


-06 38 


30.62 


114,96 


18 37 16.916 


-06 38 30 


58 


60,60 




BDS12 


4 


18 37 16.922 


-06 38 


31.57 


114,96 












BDS12 


5 


18 37 16.957 


-06 38 


30.21 


114,96 












BDS12 


6 


18 37 16.906 


-06 38 


30.48 


114,96 












BDS12 


7 


18 37 16.898 


-06 38 


30.27 


114,96 


18 37 16.897 


-06 38 30 


29 


60,60 




BDS12 


8 


18 37 16.866 


-06 38 


30.58 


114,96 












BDS12 


9 


18 37 16.875 


-06 38 


30.05 


114,96 











TABLE 2 

Integrated Flux Densities of Detected Infrared Sources 



Aperture 

Infrared Source Name 2.12 fim 8.6 (im 18.3 fim Radius 
(mjy) (Jy) (Jy) (") 



185:BDS12 1 


294 


21.9 


45.7 


1.6 


BDS12 2 


1.41 


0.12 


1.21 


0.4 


127:BDS12 1 


81.0 


7.04 


25.7 


2.6 


324:BDS12 1 




2.41 


4.67 


0.7 


BDS12 2 




0.29 


1.97 


0.5 


BDS12 3 




0.04 


<0.79 


0.4 


Combined 




2.91 


7.99 


1.5 


105:BDS12 1 


1.83 






0.2 


BDS12 2 


0.92 






1.5 


BDS12 3 


0.28 


0.28 


1.91 


0.4 


BDS12 4 


0.21 






0.2 


BDS12 5 


0.04 






0.2 


BDS12 6 


0.15 






0.1 


BDS12 7 


0.14 


0.05 


0.20 


0.1 


BDS12 8 


0.17 






0.2 


BDS12 9 


0.05 






0.1 


Combined 


6.32 


0.36 


2.02 


1.2 



Note. — All flux densities are background subtracted and derived using 
circular aperture photometry. The aperture radius given for each source 
was used for all wavelengths. "Combined" fluxes are for an aperture en- 
closing all emission in the region around the mascrs. Flux errors (lcr) are 
5% for 2.12 (Ltm, 8% for 8.6 fj,m, and 15% for 18.3 fj,m. No 2.12 ^m data 
were taken for G24. 634-00. 324. Blank entries denote no detection. Upper 
limit for detection of a point source is 2.1 mjy at 8.6 (im, and 36.7 mjy 
for 18.3 fim. 
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TABLE 3 

Auxiliary Inputs to the SED models" 



Band G23.389+00.185 G23. 657-00. 127 G24.634-00.324 G25.411+00.105 









w 


BDS12 1 

(Jy) 


BDS12 1 

(Jy) 


Combined 

(Jy) 


Combined 

(Jy) 


2MASS 


J 




1.25 


0.0007±0. 00007 


0.0005±0. 00005 


0.000767±0. 00007 


0.0002±0.0002 




H 




1.65 


0.021±0.002 


0.01±0.001 


0.005±0.0005 


0.0006±0. 00006 




K 




2.17 


0.28±0.028 


0.07±0.007 


0.064±0.006 


0.0048±0. 00048 


Spitzcr 


IRAC 


1 


3.6 


2.08±0.02 


1.2±0.02 


0.58±0.06 


0.09±0.009 




IRAC 


2 


4.5 


3.45±0.1 


2.1±0.1 


1.3±0.13 


0.4±0.04 




IRAC 


3 


5.8 


14.0± 0.4 


7.4±0.4 


3.1±0.31 


0.85±0.085 




IRAC 


4 


8.0 


_b 


7.3±1 


3.0±0.30 


1.0±0.1 


MSX 


A 




8.28 


23.47±2.4 


8.518±0.85 


4.50 


2.34 




C 




12.13 


39.00±3.9 


17.27±1.7 


6.70 


3.44 




D 




14.64 


50.99±5.1 


30.45±3.05 


8.78 


4.74 




E 




21.34 


54.71±5.5 


43.87±3.4 


13.47 


8.44 


IRAS 






12 


35.84 


16.04 


10.34 


3.54 








25 


76.44 


70.34 


25.04 


17.54 








60 


2504 


3834 


91.84 


1104 








100 


2584 


4024 


2584 


5674 


^adopted 


(kpc) 






4.47 


3.19 d 


3.0 


5.35 



Note. — All sources have IRAS and MSX detections, however if these measured flux densities arc used as upper 
limits in the SED modeling, this is denoted by a J, following those values in the table. 

a Also included in the modeling were the relevant source flux densities from Gemini observations in Table 2. 

Source is saturated in this image. 
c Assumed distance to source adopted in this paper. With the exception of G23. 657-00. 127, these arc the near 
kinematic distances for the sources based on the calculations from Rcid ct al. (2009). For SED models, the distance 
range used was this value ±10%. 

d Distance based on the trigonometric parallax (Bartkiewicz et al. 2008). 
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TABLE 4 

Derived Physical Parameters from SED Models 





^disk Rdisk 


Mi, 




Ltot 


"cavity 




(°) (AU) 


(Mq) 


(1O 3 L ) 


V 1 




G23.389+00.185:BDS12 1 








Maser ring: i=54 c 


'; _R=425AU 




Best 


18 58 


15 




15 


5 


Mode 


18 48 


14 




17 


9 


Median 


18 36 


14 




12 


6 


Range 


18-57 a -93 


13-21 




9-17 


4-25 




G23. 657-00. 127:BDS12 1 








Maser ring: i = 16 c 


'; R = 


425 AU 




Best 


18 21 


12 




6 


7 


Modc & 


18 0° 










ivieQian 


18 7 


13 




5 


a 
O 


Ron ar^ 

L V.CLIlg,'-' 


18-18 a -21 


11-16 




4-7 


5-8 


Best c 


70 179 


11 




8 


49 


Mode c 


76 179 


11 




8 


49 


Median c 


66 179 


11 




8 


49 


Range 6 


41-76 62-179 


11-12 




8-11 


^o— 00 




G24.634-00.324:Combined 








Maser ring: i = 71 c 


'; R = 


135 AU 




Best 


57 191 


9 




4 


50 


Mode d 


63 191 


9 




4 


50 


Median d 


63 217 


8 




3 


34 


Range d 


41-76 65-1880 


7-9 




3-4 


24-50 




G25.411+00.105:Combined 








Maser ring: i = 47 c 


'; R = 


550 AU 




Best 


32 12 


9 




1 


6 


Mode e 


32 


10 








Median 6 


32 157 


8 




2 


6 


Range 6 


32-57 7-782 


8-11 




1-8 


2-34 



a Means diskless. 



Results arc bimodal. Two top tens tabulated. Best overall fit 
comes from this group of results and implies a very small disk. 
Four more of these top 10 fits arc diskless. 

c Results are bimodal. Not including the best overall fit and the 
low-inclination cases. 

Results are trimodal but only showing model fits consistent 
with observations. Therefore, not showing results for those with 
disk-only fits and/or with facc-on disk inclinations. 
e One of the top ten fits was for a 1.7 Mq source and was removed. 
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Model Cross-Section 



IR bright outflow cavities 




envelope 



Scenario 1 Scenario 2 Scenario 3 




i- 0° (face-on) i~ 90° tedge-on) f~45 e 



Fig. 1. — A cross-section view of the toy model for a MYSO and 3D renderings of the model at disk/system inclinations given by the 
scenarios outlined in Section 3.2. In the cross-section view, the IR bright cavities are shown in blue and red, for what will be the blue-shifted 
and red-shifted outflow cavities, respectively, in the 3D renderings below it. The envelope and flared accretion disk are labeled. In the 
three 3D renderings of this model for the three inclination scenarios, the disk and envelope are not 100% opaque, and thus the lightening 
of the red and blue cavities reflect the effects of obscuration by the disk and envelope. In reality, the more optically thick the disk and 
envelope, the more attenuated the emission from the cavities, especially the red-shifted one. The maser rings are represented by the green 
circles, demonstrating their location with respect to the IR emission as a function of system inclination. 
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100 -10 200 -200 

ADec (mas) ADec (mas) 



G 24. 63 4-00. 3 24 



50 " 



■e 
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ADec (mas) 



200 



H 100 



3 

< 



-100 



G25.41 l-i-OO.105 
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■o - . 
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100 -100 -200 
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Fig. 2. — The four methanol maser rings in the sample with their implied disk properties, assuming the masers arise from within 
circumstellar disks. In each panel the multi-colored circles represent the individual maser spots, with colors indicating line of sight velocity 
of each spot (see Bartkiewicz et al. 2009 for these values). The dotted line is the best fitted ellipse to the maser spot distribution, with 
the cross demarcating the center of the ellipse. A toy model (as described in Figure 1) is displayed in the upper right corner of each panel, 
using the properties (from Bartkiewicz et al. 2009) given in the lower right corner. The properties given are disk inclination, position angle 
of the disk rotational symmetry axis (or the outflow axis) from north to east, and the radius of the disk given by the semi-major axis of the 
maser ring fit. The blue cone of the disk model shows the assumed blue-shifted outflow side of the MYSO, and the red cone the red-shifted 
side. The green ellipse shows the expected location of the maser ring if it is in the disk. 
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l_l i i i i i_l i i i i I i ^^J^ , , l_l , , , , I , , „ ^ _i I i , , , L 

Id"^ H.tS 1 I*.*!)' 14J5 9 |i.30 B HJJ* 1*,+B* l*4i* 1*-Sf* H.4S 1 l+.*D s 14.5BT 1 



Fig. 3. — The NIR and MIR emission registered towards G23. 389+00. 185. The top left box shows a wider field view of the 2 fim image, 
and a unsharp mask of the same image to its right. We detect a jet-like feature in the NIR that appears to be coming from source 2 (see 
appendix). The box in the top left panel is the field of view of the three panels below. The 6.7 GHz methanol masers from Figure 2 are 
represented by white crosses. The ±lcr astrometric errors in Table 1 are shown for each wavelength with the black cross in the lower left 
of each panel, along with the effective resolution (azimuthally averaged FWHM of the PSF stars) as the circle. For this source only, the 
2 fim image was sufficiently elongated that we show the 50% level of an actual star on the field in the bottom left panel. The colors on 
the bottom three panels trace the emission levels as follows: 4, 7, 20, 50, 100, 300, 900 fijy pixel -1 (at 2 fim); 2, 4, 10, 60, 150, 400, 
600 mjy pixel -1 (at 8.6 fim); and 8, 20, 60, 120, 200, 400, 600 mjy pixel -1 (at 18.3 fim). 
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Fig. 4. — The NIR and MIR emission registered towards G23. 657— 00.127. In the right panel we overlay the 2 /mi contours on the 18.3 /mi 
image. Also plotted are two lines to demonstrate how straight and flat the sides of the NIR emission are. Symbols are the same as in 
Figure 3. The colors on the bottom three panels trace the emission levels as follows: 0.1, 1, 5, 10, 20, 200, 800 (ijy pixel" 1 (at 2 /mi); 5, 
25, 50, 150, 200, 250, 300 mjy (at 8.6 /an); and 35, 75, 150, 200, 250, 300, 350 mjy pixel" 1 (at 18.3 /mi). 
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1 3 37" 22.75" 



22.70 s 22.65* 



22,50 s 22.75° 



22.70° 22.65* 



Fig. 5. — The MIR emission towards G24. 634— 00.324. The two top panels show the MIR images with their native resolution, whereas 
the bottom two panels show the same images deconvolved. Symbols are the same as in Figure 3. The colors trace the emission levels as 
follows: 0.8, 1.3, 4.7, 7.0, 14, 25, 60 mjy pixel" 1 (for 8.6 /an); 8.0, 13, 18, 23, 32, 50, 70 mjy pixel" 1 (for 18.3 £tm); 0.5, 1.3, 3.7, 10, 27, 
74, 200 mjy pixel -1 (for deconvolved 8.6 ftm); and 2.0, 4.5, 10, 22, 50, 112, 230 mjy pixel" 1 (for deconvolved 18.3 /im). 
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Fig. 6. — The NIR and MIR emission registered towards G25. 411+00. 105. The upper panel shows a wider field view of the 2 (im image 
at the natural resolution using the AO system, with the box denoting the field in the lower panels. The lower panels show the deconvolved 
images at 2, 8.6, and 18.3 (im. The 8.6 and 18.3 /an images have the 2 fim contours overlaid. Black crosses denote the masers for this source 
as given in Figure 2. All other symbols are as described in Figure 3. The colors on the bottom three panels trace the emission levels as 
follows: 0.3, 0.6, 1, 2, 4, 5, 7 ^Jy (at 2 ^m); 1, 2, 5, 10, 20, 35, 50 mjy (at 8.6 /im); and 2, 10, 25, 50, 90, 130, 170 mjy (at 18.3 fim). 
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>. ,,.rr) 



Fig. 7. — Spectral energy distributions of the continuum emission towards the four maser rings based on the data in Table 4. The filled 
circles and triangles show the input fluxes and upper limits, respectively. Red circles represent the Gemini data presented here. The black 
line shows the best fit, and the gray lines show the subsequent 9 good fits. For G23. 657-00. f 27 the results were bimodal (see Section 3.4) 
and the red fits correspond to the group of fits with i(K s fc~18 , and the blue correspond to the group of fits with idisfe ~70° . The green 
fits for G24. 634-00. 324 denote that these resultant fits are not the overall top ten best fits, but the ten best model fits after some selection 
criteria were employed (see Section 3.4). 



